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I Abstract: Absolute luminosity measurements are of general interest for colliding-beam experi- 

ments at storage rings. These measurements are necessary to determine the absolute cross-sections 
^ of reaction processes and are valuable to quantify the performance of the accelerator. Using data 

I taken in 2010, LHCb has applied two methods to determine the absolute scale of its luminosity 

Q measurements for proton-proton collisions at the LHC with a centre-of-mass energy of 7 TeV. In 

, addition to the classic "van der Meer scan" method a novel technique has been developed which 

makes use of direct imaging of the individual beams using beam-gas and beam-beam interactions. 
This beam imaging method is made possible by the high resolution of the LHCb vertex detector and 
^ the close proximity of the detector to the beams, and allows beam parameters such as positions, 

^ angles and widths to be determined. The results of the two methods have comparable precision 

CN and are in good agreement. Combining the two methods, an overal precision of 3.5% in the abso- 

O lute luminosity determination is reached. The techniques used to transport the absolute luminosity 

^ calibration to the full 2010 data-taking period are presented. 
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1. Introduction 

Absolute luminosity measurements are of general interest to coUiding-beam experiments at stor- 
age rings. Such measurements are necessary to determine the absolute cross-sections of reaction 
processes and to quantify the performance of the accelerator. The required accuracy on the value 
of the cross-section depends on both the process of interest and the precision of the theoretical pre- 
dictions. At the LHC, the required precision on the cross-section is expected to be of order 1-2%. 
This estimate is motivated by the accuracy of theoretical predictions for the production of vector 
bosons and for the two-photon production of muon pairs [1-4]. 

In a cyclical collider, such as the LHC, the average instantaneous luminosity of one pair of 
colliding bunches can be expressed as [5] 



L = NiN2f\j{vi-V2Y-- — -^^^ J pi{x,y,z,t)p2{x,y,z,t) dxdydzdt , (1.1) 

where we have introduced the revolution frequency / (11245 Hz at the LHC), the numbers of 
protons A^i and in the two bunches, the corresponding velocities vi and V2 of the particles,^ and 
the particle densities for beam 1 and beam 2, pi 2ix,y,z,t). The particle densities ai^e normalized 
such that their individual integrals over all space are unity. For highly relativistic beams colliding 
with a very small half crossing-angle a, the M0ller factor (vi — V2)^ — (vi x V2Y jc^ reduces to 
2c cos^ a ~ 2c. The integral in Eq. [O] is known as the beam overlap integral. 

Methods for absolute luminosity determination are generally classified as either direct or indi- 
rect. Indirect methods are e.g. the use of the optical theorem to make a simultaneous measurement 
of the elastic and total cross-sections [6,7], or the comparison to a process of which the absolute 
cross-section is known, either from theory or by a previous direct measurement. Direct measure- 
ments make use of Eq. |0] and employ several strategies to measure the various parameters in the 
equation. 

The analysis described in this paper relies on two direct methods to determine the absolute 
luminosity calibration: the "van der Meer scan" method (VDM) [8,9] and the "beam-gas imaging" 
method (BGI) [10]. The BGI method is based on reconstructing beam-gas interaction vertices to 
measure the beam angles, positions and shapes. It was applied for the first time in LHCb (see 
Refs. [1 1-13]) using the first LHC data collected at the end of 2009 at ^ = 900 GeV. The BGI 
method relies on the high precision of the measurement of interaction vertices obtained with the 
LHCb vertex detector. The VDM method exploits the ability to move the beams in both transverse 
coordinates with high precision and to thus scan the colliding beams with respect to each other. 
This method is also being used by other LHC experiments [14]. The method was first applied 
at the CERN ISR [8]. Recently it was demonstrated that additional information can be extracted 
when the two beams probe each other such as during a VDM scan, allowing the individual beam 
profiles to be determined by using vertex measurements of pp interactions in beam-beam collisions 
(beam-beam imaging) [15]. 

In principle, beam profiles can also be obtained by scanning wires across the beams [16] or 
by infening the beam properties by theoretical calculation from the beam optics. Both methods 
lack precision, however, as they both rely on detailed knowledge of the beam optics. The wire-scan 

'in the approximation of zero emittance the velocities are the same within one bunch. 
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method is limited by the achievable proximity of the wire to the interaction region which introduces 
the dependence on the beam optics model. 

The LHC operated with a pp centre-of-mass energy of 7 TeV (3.5 TeV per beam). Typical 
values observed for the transverse beam sizes are close to 50 (xm and 55 mm for the bunch length. 
The half-crossing angle was typically 0.2 mrad. 

Data taken with the LHCb detector, located at interaction point (IP) 8, are used in conjunction 
with data from the LHC beam instrumentation. The measurements obtained with the VDM and 
BGI methods are found to be consistent, and an average is made for the final result. The limiting 
systematics in both measurements come from the knowledge of the bunch populations A^i and N2- 
All other sources of systematics are specific to the analysis method. Therefore, the comparison of 
both methods provides an important cross check of the results. The beam-beam imaging method is 
applied to the data taken during the VDM scan as an overall cross check of the absolute luminosity 
measurement. 

Since the absolute calibration can only be performed during specific running periods, a relative 
normalization method is needed to transport the results of the absolute calibration of the luminos- 
ity to the complete data-taking period. To this end we defined a class of visible interactions. The 
cross-section for these interactions is determined using the measurements of the absolute luminos- 
ity during specific data-taking periods. Once this visible cross-section is determined, the integrated 
luminosity for a period of data-taking is obtained by accumulating the count rate of the corre- 
sponding visible interactions over this period. Thus, the calibration of the absolute luminosity is 
translated into a determination of a well defined visible cross-section. 

In the present paper we first describe briefly the LHCb detector in Sect. @, and in particular 
those aspects relevant to the analysis presented here. In Sect. the methods used for the relative 
normalization technique are given. The determination of the number of protons in the LHC bunches 
is detailed in Sect. 3. The two methods which are used to determine the absolute scale are described 
in Sect. ^ and ^ respectively. The cross checks made with the beam-beam imaging method are 
shown in Sect. 0. Finally, the results are combined in Sect. |[ 

2. The LHCb detector 

The LHCb detector is a magnetic dipole spectrometer with a polar angular coverage of approxi- 
mately 10 to 300 mrad in the horizontal (bending) plane, and 10 to 250 mrad in the vertical plane. 
It is described in detail elsewhere [17]. A right-handed coordinate system is defined with its origin 
at the nominal pp interaction point, the z axis along the average nominal beam line and pointing 
towards the magnet, and the y axis pointing upwards. Beam 1 (beam 2) travels in the direction of 
positive (negative) z- 

The apparatus contains tracking detectors, ring-imaging Cherenkov detectors, calorimeters, 
and a muon system. The tracking system comprises the vertex locator (VELO) suiTounding the 
pp interaction region, a tracking station upstream of the dipole magnet and three tracking stations 
located downstream of the magnet. Particles traversing the spectrometer experience a bending-field 
integral of around 4 Tm. 

The VELO plays an essential role in the application of the beam-gas imaging method at LHCb. 
It consists of two retractable halves, each having 21 modules of radial and azimuthal silicon-strip 
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beam1 



Figure 1. A sketch of the VELO, including the two Pile-Up stations on the left. The VELO sensors are 
drawn as double lines while the PU sensors are indicated with single lines. The thick arrows indicate the 
direction of the LHC beams (beam 1 going from left to right), while the thin ones show example directions 
of flight of the products of the beam-gas and beam-beam interactions. 

sensors in a half-circle shape, see Fig. [I]. Two additional stations {Pile-Up System, PU) upstream 
of the VELO tracking stations are mainly used in the hardware trigger. The VELO has a large 
acceptance for beam-beam interactions owing to its many layers of silicon sensors and their close 
proximity to the beam line. During nominal operation, the distance between sensor and beam is 
only 8 mm. During injection and beam adjustments, the two VELO halves are moved apart in a 
retracted position away from the beams. They are brought to their nominal position close to the 
beams during stable beam periods only. 

The LHCb trigger system consists of two separate levels: a hardware trigger (LO), which is 
implemented in custom electronics, and a software High Level Trigger (HLT), which is executed 
on a farm of commercial processors. The LO trigger system is designed to run at 1 MHz and uses 
information from the Pile-Up sensors of the VELO, the calorimeters and the muon system. They 
send information to the LO decision unit (LODU) where selection algorithms are run synchi^onously 
with the 40 MHz LHC bunch-crossing signal. For every nominal bunch-crossing slot {i.e. each 25 
ns) the LODU sends decisions to the LHCb readout supervisor. The full event information of all 
sub-detectors is available to the HLT algorithms. 

A trigger strategy is adopted to select pp inelastic interactions and collisions of the beam with 
the residual gas in the vacuum chamber. Events are collected for the four bunch-crossing types: 
two colliding bunches (bb), one beam 1 bunch with no beam 2 bunch (be), one beam 2 bunch with 
no beam 1 bunch (eb) and nominally empty bunch slots (ee). Here "b" stands for "beam" and "e" 
stands for "empty". The first two categories of crossings, which produce particles in the forward 
direction, are triggered using calorimeter information. An additional PU veto is applied for be 
crossings. Crossings of the type eb, which produce particles in the backward direction, are triggered 
by demanding a minimal hit multiplicity in the PU, and vetoed by calorimeter activity. The trigger 
for ee crossings is defined as the logical OR of the conditions used for the be and eb crossings 
in order to be sensitive to background from both beams. During VDM scans specialized trigger 
conditions are defined which optimize the data taking for these measurements (see Sect. |57T1 ). 

The precise reconstruction of interaction vertices ("primary vertices", PV) is an essential in- 
gredient in the analysis described in this paper. The initial estimate of the PV position is based on 
an iterative clustering of tracks ("seeding"). Only tracks with hits in the VELO are considered. For 
each track the distance of closest approach (DOCA) with all other tracks is calculated and tracks 
are clustered into a seed if their DOCA is less than 1 mm. The position of the seed is then obtained 
using an iterative procedure. The point of closest approach between all track pairs is calculated and 
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its coordinates are used to discard outliers and to determine the weighted average position. The 
final PV coordinates are determined by iteratively improving the seed position with an adaptive, 
weighted, least-squares fit. In each iteration a new PV position is evaluated. Participating tracks are 
extrapolated to the z coordinate of the PV and assigned weights depending on their impact parame- 
ter with respect to the PV. The procedure is repeated for all seeds, excluding tracks from previously 
reconstructed primary vertices, retaining only PVs with at least five tracks. For this analysis only 
PVs with a larger number of tracks are used since they have better resolution. For the study of 
beam-gas interactions only PVs with at least ten tracks are used and at least 25 tracks are required 
for the study of pp interactions. 

3. Relative normalization method 

The absolute luminosity is obtained only for short periods of data-taking. To be able to perform 
cross-section measurements on any selected data sample, the relative luminosity must be measured 
consistently during the full period of data taking. The systematic relative normalization of all data- 
taking periods requires specific procedures to be appUed in the trigger, data-acquisition, processing 
and final analysis. The basic principle is to acquire luminosity data together with the physics data 
and to store it in the same files as the physics event data. During further processing of the physics 
data the relevant luminosity data is kept together in the same storage entity. In this way, it remains 
possible to select only part of the full data-set for analysis and still keep the capability to determine 
the corresponding integrated luminosity. 

The luminosity is proportional to the average number of visible proton-proton interactions in 
a beam-beam crossing, jx^-^^. The subscript "vis" is used to indicate that this holds for an arbitrary 
definition of the visible cross-section. Any stable interaction rate can be used as relative luminosity 
monitor. For a given period of data-taking, the integrated interaction rate can be used to determine 
the integrated luminosity if the cross-section for these visible interactions is known. The determi- 
nation of the cross-section corresponding to these visible interactions is achieved by cahbrating the 
absolute luminosity during specific periods and simultaneously counting the visible interactions. 

Triggers which initiate the full readout of the LHCb detector are created for random beam 
crossings. These are called "luminosity triggers". During normal physics data-taking, the overall 
rate is chosen to be 997 Hz, with 70% assigned to bb, 15% to be, 10% to eb and the remaining 
5% to ee crossings. The events taken for crossing types other than bb are used for background 
subtraction and beam monitoring. After a processing step in the HLT a small number of "luminosity 
counters" are stored for each of these random luminosity triggers. The set of luminosity counters 
comprise the number of vertices and tracks in the VELO, the number of hits in the PU and in the 
scintillator pad detector (SPD) in front of calorimeters, and the transverse energy deposition in the 
calorimeters. Some of these counters are directly obtained from the LO, others are the result of 
partial event-reconstruction in the HLT. 

During the final analysis stage the event data and luminosity data are available on the same 
files. The luminosity counters are summed (when necessary after time-dependent calibration) and 
an absolute calibration factor is applied to obtain the absolute integrated luminosity. The absolute 
calibration factor is universal and is the result of the luminosity calibration procedure described in 
this paper. 
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The relative luminosity can be determined by summing the values of any counter which is 
linear with the instantaneous luminosity. Alternatively, one may determine the relative luminos- 
ity from the fraction of "empty" or invisible events in bb crossings which we denote by P^. An 
invisible event is defined by applying a counter-specific threshold below which it is considered 
that no pp interaction was seen in the corresponding bunch crossing. Since the number of events 
per bunch crossing follows a Poisson distribution with mean value proportional to the luminosity, 
the luminosity is proportional to — ItiPq. This "zero count" method is both robust and easy to im- 
plement [18]. In the absence of backgrounds, the average number of visible pp interactions per 
crossing can be obtained from the fraction of empty bb crossings by /i^jj, = — InP^^. Backgrounds 
are subtracted using 

= - (inP^^ - InP^' - InPt + InPo^^) , (3.1) 

where PQ{i = bb, ee, be, eb) are the probabilities to find an empty event in a bunch-crossing slot for 
the four different bunch-crossing types. The Pq^ contribution is added because it is also contained 
in the Pq^ and Pq'" terms. The purpose of the background subtraction, Eq. ^j] , is to correct the 
count-rate in the bb crossings for the detector response which is due to beam-gas interactions and 
detector noise. In principle, the noise background is measured during ee crossings. In the presence 
of parasitic beam protons in ee bunch positions, as will be discussed below, it is not correct to 
evaluate the noise from P^^. In addition, the detector signals are not fully confined within one 
25 ns bunch-crossing slot. The empty (ee) bunch-crossing slots immediately following a bb, be or 
eb crossing slot contain detector signals from interactions occurring in the preceding slot ("spill- 
over"). The spill-over background is not present in the bb, be and eb crossings. Therefore, since 
the detector noise for the selected counters is small (< 10^^ relative to the typical values measured 
during bb crossings) the term InP^^ in Eq. ^7T] is neglected. Equation ^7T] assumes that the proton 
populations in the be and eb crossings are the same as in the bb crossings. With a population spread 
of typically 10% and a beam-gas background fraction < 10^^ compai^ed to the pp interactions the 
effect of the spread is negligible, and is not taken into account. 

The results of the zero-count method based on the number of hits in the PU and on the number 
of tracks in the VELO are found to be the most stable ones. An empty event is defined to have 
< 2 hits when the PU is considered or < 2 tracks when the VELO is considered. A VELO track 
is defined by at least three hits on a straight line in the radial strips of the silicon detectors of the 
VELO. The number of tracks reconstructed in the VELO is chosen as the most stable counter. In 
the following we will use the notation a^j^, (= CJygLo) for the visible cross-section measured using 
this method, except when explicitly stated otherwise. Modifications and alignment variations of the 
VELO also have negligible impact on the method, since the efficiency for reconstructing at least 
two tracks in an inelastic event is very stable against detector effects. Therefore, the systematics 
associated with this choice of threshold is negligible. The stability of the counter is demonstrated 
in Fig. which shows the ratio of the relative luminosities determined with the zero-count method 
from the multiplicity of hits in the PU and from the number of VELO tracks. Apart from a few 
threshold updates in the PU configuration, the PU was also stable throughout LHCb 2010 running, 
and it was used as a cross check. Figure covers the whole period of LHCb operation in 2010, 
with both low and high number of interactions per crossing. Similar cross checks have been made 
with the counter based on the number of reconstructed vertices. These three counters have different 
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Figure 2. Ratio between /i^j^ values obtained with the zero-count method using the number of hits in the 
PU and the track count in the VELO versus /^velo- deviation from unity is due to the difference in 
acceptance. The left (right) panel uses runs from the beginning (end) of the 2010 running period with lower 
(higher) values of Mvelo- horizontal lines indicate a ±1% variation. 

systematics, and by comparing their ratio as a function of time and instantaneous luminosity we 
conclude that the relative luminosity measurement has a systematic error of 0.5%. 

The number of protons, beam sizes and transverse offsets at the interaction point vary across 
bunches. Thus, the /x^jj, value varies across bb bunch crossings. The spread in /x^jj, is about 10% 
of the mean value for typical runs. Due to the non-linearity of the logarithm function one first 
needs to compute fj.^^^ values for different bunch crossings and then to take the average. However, 
for short time intervals the statistics are insufficient to distinguish between bunch-crossing IDs, 
while one cannot assume fi^^^ to be constant when the intervals are too long due to e.g. loss of 
bunch population and emittance growth. If the spread in instantaneous /x^jj, is known, the effect of 
neglecting it in calculating an average value of /x^jj, can be estimated. The difference between the 
naively computed {jl^^^ value and the true one is then 

= -Hn) - (-(ta^o)) = (MtIv)) - (3-2) 

where the average is taken over all beam-beam crossing slots /. Therefore, the biased /i^jj value 
can be calculated over short time intervals and a correction for the spread of /x^jj, can in principle 
be applied by computing Pq/{Pq) over long time intervals. At the present level of accuracy, this 
correction is not required.^ The effect is only weakly dependent on the luminosity counter used. 

4. Bunch population measurements 

To measure the number of particles in the LHC beams two types of beam current transformers are 
installed in each ring [19]. One type, the DCCT (DC Current Transformer), measures the total 
current of the beams. The other type, the FBCT (Fast Beam Current Transformer), is gated with 

■^The relative luminosity increases by 0.5% when the correction is applied. 
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Figure 3. Ratio between jU^j^, values obtained with the zero-count method using the number of hits in the PU 
and the track count in the VELO as a function of time in seconds relative to the first run of LHCb in 2010. 
The period spans the full 2010 data taking period (about half a year). The dashed lines show the average 
value of the starting and ending periods (the first and last 25 runs, respectively) and differ by w 1%. The 
changes in the average values between the three main groups (t < 0.4 x 10^ s, 0.4 x 10^ < t < 1.2 x 10^ s, 
f > 1.2 X 10^ s) coincide with known maintenance changes to the PU system. The upward excursion near 
1.05 X 10^ s is due to background introduced by parasitic collisions located at 37.5 m from the nominal IP 
present in the bunch filling scheme used for these fills to which the two counters have different sensitivity. 
The downward excursion near 0.25 x 10^ s is due to known hardware failures in the PU (recovered after 
maintenance). The statistical errors are smaller than the symbol size of the data points. 



25 ns intervals and is used to measure the relative charges of the individual bunches. The DCCT is 
absolutely calibrated, thus is used to constrain the total number of particles, while the FBCT defines 
the relative bunch populations. The procedure is described in detail in Ref. [20]. All devices have 
two independent readout systems. For the DCCT both systems provide reliable information and 
their average is used in the analysis, while for the FBCT one of the two systems is dedicated to 
tests and cannot be used. 

The absolute calibration of the DCCT is determined using a high-precision current source. 
At low intensity (early data) the noise in the DCCT readings is relatively important, while at the 
higher intensities of the data taken in October 2010 this effect is negligible. The noise level and 
its variation is determined by interpolating the average DCCT readings over long periods of time 
without beam before and after the relevant fills. 

In addition to the absolute calibration of the DCCTs, a deviation from the proportionality of 
the FBCT readings to the individual bunch charges is a potential source of systematic uncertainty. 
The FBCT charge offsets are cross checked using the ATLAS BPTX (timing) system [21]. This 
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comparison shows small discrepancies between their offsets. These deviations are used as an esti- 
mate of the uncertainties. Since the FBCT equipment is readjusted at regular intervals, the offsets 
can vary on a fill-by-fill basis. Following the discussion in Ref. [20], an estimate of 2.9% is used 
for the uncertainty of an individual bunch population product of a colliding bunch pair. Owing to 
the DCCT constraint on the total beam current, the overal uncertainty is reduced when averaging 
results of different bunch pairs within a single fill. As will be discussed in Sect. 0, for the analysis 
of the VDM data a method can be used which only needs the assumption of the linearity of the 
FBCT response. 

The LHC radio frequency (RF) system operates at 400 MHz, compared to the nominal 40 MHz 
bunch frequency. If protons circulate in the ring outside the nominal RF buckets, the readings of 
the DCCT need to be corrected before they are used to normalize the sum of the FBCT signals. We 
define "satelUte" bunches as charges in neighbouring RF buckets compared to the nominal bucket. 
Satellite bunches inside the nominally filled bunch slots can be detected by the LHC experiments 
when there is no (or a very small) crossing angle between the two beams. The satellites would be 
observed as interactions displaced by a multiple of 37.5 cm from the nominal intersection point. 
For a part of the 2010 run the ATLAS and CMS experiments were operating with zero crossing 
angle and displaced interactions were indeed observed [20]. 

The "ghost charge" is defined as the charge outside the nominally filled bunch slots. The 
rates of beam-gas events produced by "ghost" and nominal protons are measured using the beam- 
gas trigger. The ghost fraction is determined by comparing the number of beam-gas interactions 
during ee crossings with the numbers observed in be and eb crossings. The timing of the LHCb 
trigger is optimized for interactions in the nominal RF buckets. The trigger efficiency depends on 
the time of the interaction with respect to the phase of the clock (modulo 25 ns). A measurement 
of the trigger efficiency was performed by shifting the clock which is usually synchronized with 
the LHC bunch-crossing time by 5, 10 and 12.5 ns and by comparing the total beam-gas rates 
in the nominal crossings. From these data the average efficiency for ghost charge is obtained to 
be Eaverage = 0.86 it 0.14 (0.84 it 0.16) for beam 1 (beam 2). The ghost charge is measured for 
each fill during which an absolute luminosity measurement is performed and is typically 1% of 
the total beam charge or less. The contribution of "ghost" protons to the total LHC beam current 
is subtracted from the DCCT value before the sum of the FBCT bunch populations is constrained 
by the DCCT measurement of the total current. The uncertainty assigned to the subtraction of 
ghost charge varies per fill and is due to the trigger efficiency uncertainty and the limited statistical 
accuracy. These two error components are of comparable size. 

5. The van der Meer scan (VDM) method 

The beam position scanning method, invented by van der Meer, provides a direct determination of 
an effective cross-section a^-j, by measuring the corresponding counting rate as a function of the 
position offsets of two colliding beams [8]. At the ISR only vertical displacements were needed 
owing to the crossing angle between the beams in the horizontal plane and to the fact that the beams 
were not bunched. For the LHC the beams have to be scanned in both transverse directions due to 
the fact that the beams are bunched [9]. The cross-section a^^^ can be measured for two colliding 
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Table 1. Parameters of LHCb van der Meer scans. A^i 2 is the typical number of protons per bunch, j3* 
characterizes the beam optics near the IP, «tot (wcoii) is the total number of (colliding) bunches per beam, 
/^vTs" '■^^ average number of visible interactions per crossing at the beam positions with maximal rate. 
^2 is the decay time of the product of the bunch populations and Tj, is the decay time of the luminosity. 

25 Apr 15 Oct 



LHC fill number 
Ni,2 (10^° protons) 
r (m) 

"coll/«tot 
..max 

Trigger 



T^NiNi (h) 

TL(h) 



1059 
1 
2 

1/2 
0.03 
minimum bias 



950 
30 



1422 
7-8 
3.5 
12/16 
1 

22.5 kHz random 
430 Hz minimum bias 
beam-gas 
700 
46 



bunches using the equation [15] 

^ Mis (Ax, Ay J JA,- Jn^.^{A,^,Ay)dAy 

MA^2Mvis(Avo,A,o)cosa ' ^ " ^ 

where /x^;j,(A_r, Ay) is the average number of interactions per crossing at offset {Ax, Ay) correspond- 
ing to the cross-section a^j^,. The interaction rates R{Ax,Ay) are related to iJ.^^^{Ax,Ay) by the rev- 
olution frequency, R{Ax,Ay) = /Mvis(Ax,Av). These rates are measured at offsets A^^ and Ay with 
respect to their nominal positions at offsets {Axg,Ayg). The scans consist of creating offsets Av and 
Ay such that practically the full profiles of the beams are explored. The measured rate integrated 
over the displacements gives the cross-section. 

The main assumption is that the density distributions in the orthogonal coordinates x and y can 
be factorized. In that case, two scans are sufficient to obtain the cross-section: one along a constant 
j-displacement Ay^ and one along a constant ;c-displacement A^^. It is also assumed that effects due 
to bunch evolution during the scans (shape distortions or transverse kicks due to beam-beam effects, 
emittance growth, bunch current decay), effects due to the tails of the bunch density distribution 
in the transverse plane and effects of the absolute length scale calibration against magnet current 
trims are either negligible or can be corrected for. 

5.1 Experimental conditions during the van der Meer scan 

VDM scans were performed in LHCb during dedicated LHC fills at the beginning and at the end 
of the 2010 running period, one in April and one in October. The characteristics of the beams are 
summarized in Table [I]. In both fills there is one scan where both beams moved symmetrically 
and one scan where only one beam moved at a time. Precise beam positions are calculated from 
the LHC magnet currents and cross checked with vertex measurements using the LHCb VELO, as 
described below. 
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In April the maximal beam movement of ±3 a was achieved only in the first scan, as in the 
second only the first beam was allowed to move.^ During the second October scan, both beams 
moved one after the other, covering the whole separation range of ^ 6a to both sides. However, 
the beam steering procedure was such that in the middle of the scan the first beam jumped to 
an opposite end point and then returned, so that the beam movement was not continuous. This 
potentially increases hysteresis effects in the LHC magnets. In addition, the second scan in October 
had half the data points, so it was used only as a cross check to estimate systematic errors. 

During the April scans the event rate was low and it was possible to record all events containing 
visible interactions. A loose minimum bias trigger was used with minimal requirements on the 
number of SPD hits (> 3) and the transverse energy deposition in the calorimeters (> 240 MeV). In 
October the bunch populations were higher by a factor ~ 7.5, therefore, in spite of slightly broader 
beams (the optics defined a j8* value of 3.5 m instead of 2 m in April), the rate per colliding bunch 
pair was higher by a factor of ~ 30. There were twelve colliding bunch pairs instead of one in 
April. Therefore, a selective trigger was used composed of the logical OR of three independent 
criteria. The first decision accepted random bunch crossings at 22.5 kHz (20 kHz were devoted 
to the twelve crossings with collisions, 2 kHz to the crossings where only one of two beams was 
present, and 0.5 kHz to the empty crossings). The second decision used the same loose minimum 
bias trigger as the one used in April but its rate was limited to 130 Hz. The third decision collected 
events for the beam-gas analysis. 

For both the April and October data the systematic error is dominated by uncertainties in the 
bunch populations. In April this uncertainty is higher (5.6%) due to a larger contribution from 
the offset uncertainty at lower bunch populations [20] . In October the measurement of the bunch 
populations was more precise, but its uncertainty (2.7%) is still dominant in the cross-section de- 
termination [22]. Since the dominant uncertainties are systematic and correlated between the two 
scans, we use the less precise April scan only as a cross check. The scans give consistent results, 
and in the following we concentrate on the scan taken in October which gives about a factor two 
better overall precision in the measurement. 

The LHC filling scheme was chosen in such a way that all bunches collided only in one exper- 
iment (except for ATLAS and CMS where the bunches are always shared), namely twelve bunch 
pairs in LHCb, three in ATLAS/CMS and one in ALICE. The populations of the bunches colliding 
in LHCb changed during the two LHCb scans by less than 0.1%. Therefore, the rates are not nor- 
malized by the bunch population product A^i of each colliding bunch pair at every scan point, 
but instead only the average of the product over the scan duration is used. This is done to avoid the 
noise associated with the M,2 measurement. The averaged bunch populations are given in Table ^ 
The same procedure is applied for the April scan, when the decay time of Ni was longer, 950 
instead of 700 hours in October. 

In addition to the bunch population changes, the luminosity stability may be limited by the 
changes in the bunch profiles, e.g. by emittance growth. The luminosity stability is checked several 
times during the scans when the beams were brought back to their nominal position. The average 
number of interactions per crossing is shown in Fig. ^ for the October scan. The luminosity decay 
time is measured to be 46 hours (30 hours in April). This corresponds to a 0.7% luminosity drop 

^We refer here to la as the average of the approximate widths of the beams. 
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Table 2. Bunch populations (in 10 particles) averaged over the two scan periods in October separately. 
The bottom line is the DCCT measurement, all other values are given by the FBCT The first 12 rows are the 
measurements in bunch crossings (BX) with collisions at LHCb, and the last two lines are the sums over all 
16 bunches. 

Scan 1 Scan 2 



BX 


Ni 


N2 


Ni 


A^2 


2027 


8.425 


7.954 


8.421 


7.951 


2077 


7.949 


7.959 


7.944 


7.957 


2127 


7.457 


7.563 


7.452 


7.561 


2177 


6.589 


7.024 


6.584 


7.021 


2237 


7.315 


8.257 


7.311 


8.255 


2287 


7.451 


7.280 


7.446 


7.278 


2337 


7.016 


7.219 


7.012 


7.217 


2387 


7.803 


6.808 


7.798 


6.805 


2447 


7.585 


7.744 


7.580 


7.742 


2497 


7.878 


7.747 


7.874 


7.745 


2547 


6.960 


6.244 


6.955 


6.243 


2597 


7.476 


7.411 


7.472 


7.409 


All, FBCT 


120.32 


119.07 


120.18 


118.99 


DCCT 


120.26 


119.08 


120.10 


118.98 



during the first, longer, scan along either Av or Ay (0.9% in April). The scan points have been taken 
from lower to higher A^, Ay values, therefore, the luminosity drop effectively enhances the left part 
of the integral and reduces its right part, so that the net effect cancels to first order since the curve 
is symmetric. The count rate /?(Ayg, Ayg) at the nominal position entering Eq. ^TT], is measured in 
the beginning, in the middle and at the end of every scan, so that the luminosity drop also cancels 
to first order. Therefore, the systematic error due to the luminosity drop is much less than 0.7% and 
is neglected. 

The widths of the profiles of the luminous region did not change within the statistical uncer- 
tainties when the beams were brought to their nominal positions during the first and the second 
scans in Aj^ and Ay. In addition, the width of the profiles measured in the two VDM scans did 
not change. These facts also indicate that the effect of the emittance growth on the cross-section 
measurement is negligible. 

5.2 Cross-section determination 

In accordance with the definition of the most stable relative luminosity counter, a visible event is 
defined as a pp interaction with at least two VELO tracks. The twelve colliding bunch pairs of 
the VDM scan in October are analysed individually. The dependence on the separation Ajc and Ay 
of jUyjj, summed over all bunches is shown in Fig. ^ Two scans are overlaid, the second is taken 
at the same values of Ax and Ay but with twice as large a step size and different absolute beam 
positions. One can see that the Ay curves are not well reproduced in the two scans. The reason 
for this apparent non-reproducibility is not understood. It may be attributed to hysteresis effects or 
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Figure 4. Evolution of the average number of interactions per crossing at the nominal beam position during 
the October scans. In the first (second) scan the parameters at the nominal beam position were measured 
three (four) times both during the scan and the A, scan. The straight line is a fit to the data. The luminosity 
decay time is 46 hours. 

Table 3. Mean and RMS of the VDM count-rate profiles summed over the twelve colliding bunch pairs 
obtained from data in the two October scans (scan 1 and scan 2). The statistical errors are 0.05 i.im in the 
mean position and 0.04 |j.m in the RMS. 



imperfections in the description of the optics.^ 

The mean and RMS values of the VDM count-rate profiles shown in Fig. ^ are listed in Table ^ 
Single Gaussian fits to the individual bunch profiles return values between 2.7 and 4.3 per degree 
of freedom. Double Gaussian fits provide a much better description of the data and are therefore 
used in the analysis. The single Gaussian fits give cross-section values typically 1.5 to 2% larger 
than the ones obtained with a double Gaussian. It is found that the fit errors can be reduced by 
approximately a factor two if the fits to the and Ay curves are performed simultaneously and 
the value measured at the nominal point /x^jg ( A^q ,AyQ ) is constrained to be the same in both scans. 
The first fit parameter is chosen to be / /i^j^JAj / /i^jj,t/Ay//i^jj,(A;rp, Ay^), so that a correlation of 
both integrals and the value at the nominal point is correctly taken into account in the resulting 
fit error. Other fit parameters are: the two integrals along Ay and Ay, and 0\, Aa and a common 
central position of the Gaussian function for the A;^ and similarly for the A^, curves. Here G\ and 



Scan Ay scan Ay scan 



Mean (|J.m) 



RMS (^m) 



1 13 3T 

2 2.8 9.2 

1 80.6 80.8 

2 80.5 80.7 



''imperfections in the description of the optics can manifest themselves as second order effects in the translation of 
magnet settings into beam positions or beam angles. 




y separation {\im) 

Figure 5. Number of interactions per crossing summed over the twelve colliding bunches versus the separa- 
tions A, (top), A,, (bottom) in October. The first (second) scan is represented by the dark/blue (shaded/red) 
points and the solid (dashed) lines. The spread of the mean values and widths of the distributions obtained 
individually for each colliding pair are small compared to the widths of the VDM profiles, so that the sum 
gives a good illustration of the shape. The curves represent the single Gaussian fits to the data points de- 
scribed in the text. 



02 = Y of + Aa^ are the two Gaussian widths of the fit function. The relative normalization of 
the two Gaussian components and the value at the nominal point are derived from the nine fit 
parameters listed above. The value per degree of freedom of the fit is between 0.7 and 1.8 for 
all bunch pairs. 

The product of bunch populations Ni N2 of the twelve colliding bunches have an RMS spread 
of 12%. The analysis of the individual bunch pairs gives cross-sections consistent within statistical 
errors, which typically have values of 0.29% in the first scan. The sensitivity of the method is high 
enough that it is possible to calibrate the relative bunch populations N\ 2/ L]=i ^1 2 measured with 
the FBCT system by assuming a linear response. Here / runs over the twelve bunches colliding 
in LHCb and j over all 16 bunches circulating in the machine. By comparing the FBCT with 
the ATLAS BPTX measurements it is observed that both may have a non-zero offset [20,22]. A 
discrete function s[^^ is fitted to the twelve measurements a^j^, using three free parameters: the 
common cross-section a^^^ and the two FBCT offsets for the two beams A^J'j 



n 

b=l,2 



(5.2) 
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Table 4. Results for the visible cross-section fitted over the twelve bunches colliding in LHCb for the 
October VDM data together with the results of the April scans. N^2 ^6 the FBCT or BPTX offsets in units 
of 10^" particles. They should be subtracted from the values measured for individual bunches. The first (last) 
two columns give the results for the first and the second scan using the FBCT (BPTX) to measure the relative 
bunch populations. The cross-section from the first scan obtained with the FBCT bunch populations with 
offsets determined by the fit is used as final VDM luminosity calibration. The results of the April scans are 
reported on the last row. Since there is only one colliding bunch pair, no fit to the FBCT offsets is possible. 



October data 





FBCT 


ATLAS BPTX 




Scan 1 


Scan 2 


Scan 1 Scan 2 




with fitted offsets 


with fitted offsets 




58.73 ± 0.05 


57.50 ± 0.07 


58.62 ± 0.05 57.45 ± 0.07 




0.40 ±0.10 


0.29 ±0.15 


-0.10 ±0.12 -0.23 ±0.17 




-0.02 ±0.10 


0.23 ±0.13 


-0.63 ±0.12 -0.34 ±0.15 




5.8/9 


7.6/9 


6.9/9 7.3/9 




with offsets fixed at zero 


with offsets fixed at zero 


f^vis (mb) 


58.73 ± 0.05 


57.50 ± 0.07 


58.63 ± 0.05 57.46 ± 0.07 




23.5/11 


21.9/11 


66.5/ 11 23.5/ 11 


Apri 


1 data 




Scan 1 


Scan 2 




f^vis (mb) 


59.6±0.5 


57.0±0.5 





which corrects the relative populations ^ for the FBCT offsets N\.^ and takes into account that 
the total beam intensities measured with the DCCT constrain the sums of all bunch populations 
obtained from the FBCT values. The sum over all bunches of the quantites Ljli^ii' 
normalized to the DCCT value prior to the fit and the fit using Eq. evaluates the correction due 
to the FBCT offsets alone. The results of this fit are shown in Fig. ^ where the data points a^jj, are 
drawn without offset correction and the lines represent the fit function of Eq. The use of two 
offsets improves the description of the points compared to the uncorrected simple fit. The per 
degree of freedom and other relevant fit results are summarized in Table ^ In addition, the table 
also shows results for the case where the ATLAS BPTX is used instead of the FBCT system. 

One can see that the offset errors in the first scan are (0.10-0.12) X 10^°, or 1.5% relative 
to the average bunch population (^,2) = 7.5 x lO"'. The sensitivity of the method, therefore, is 
very high, in spite of the fact that the RMS spread of the bunch population products A/^i is 12%. 
The quoted errors are only statistical. For the first scan, the relative cross-section error is 0.09%. 
Since the fits return good x'- values, the bunch-crossing dependent systematic uncertainties (such 
as emittance growth and bunch population product drop) are expected to be lower or comparable. 
An indication of the level of the systematic eiTors is given by the difference of about two standard 
deviations found for 'N^^ between the two scans. All principal sources of systematic eiTors which 
will be discussed below (DCCT scale uncertainty, hysteresis, and ghost charges) cancel when com- 
paring bunches within a single scan. 

In spite of the good agreement between the bunches within the same scan, there is an overall 
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Figure 6. Cross-sections without correction for the FBCT offset for the twelve bunches of the October VDM 
fill (data points). The lines indicate the results of the fit as discussed in the text. The upper (lower) set of 
points is obtained in the first (second) scan. 



2.1% discrepancy between the scans. The reason is not understood, and may be attributed to a 
potential hysteresis effect or similar effects resuhing in uncontrollable shifts of the beam as a whole. 
The results of the first scan with the FBCT offsets determined by the fit are taken as the final VDM 
luminosity determination (see Sect. The 2.1% uncertainty estimated from the discrepancy is 
the second largest systematic error in the cross-section measurement after the uncertainties in the 
bunch populations. In the April data the situation is similar: the discrepancy between the cross- 
sections obtained from the two scans is (4.4 it 1 .2)%, the results may be found in Table |[ Since the 
April measurement is performed using corrected trigger rates proportional to the luminosity instead 
of VELO tracks, the results have been corrected for the difference in acceptances. The correction 
factor is determined by studying random triggers and is CTygLo/^^'Aprii trigger — 1-066, where c^ygLo 
is the usual definition of a„:„. 



5.3 Systematic errors 

5.3.1 Reproducibility of the luminosity at the nominal beam positions 

Figure ^ shows the evolution of the luminosity as a function of time for the periods where the 
beams were at their nominal positions during the VDM scan. One expects a behaviour which 
follows the loss of beam particles and the emittance growth. Since these effects occur at large 
time-scales compared to the duration of the scan, the dependence on these known effects can be 
approximated by a linear evolution. As shown in Fig. 0, the luminosity did not always return to 
the expected value when the beams returned to their nominal positions. The ;^^/ndf with respect 
to the fitted straight line is too large (40/12), thus, the non-reproducibility cannot be attributed 
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fully to statistical fluctuations and another systematic effect is present. The origin of this effect is 
not understood but it may be similar to the one which causes the non-reproducibility of the beam 
positions observed in the shift of the two scan curves. Therefore, a systematic error of 0.4% is 
assigned to the absolute scale of the jJ.^^^ measurement to take this observation into account. The 
systematic error is estimated as the amount which should be added in quadrature to the statistical 
error of 0.25% to produce a ;^^/ndf equal to one. Since the absolute scale of the jJ.^^^ measurement 
enters the cross-section linearly (Eq. |57T1 ), the same systematic en^or of 0.4% is assigned to the 
cross-section measurement. 

5.3.2 Length scale calibration 

The beam sepai^ation values and Ay ai^e calculated from the LHC magnet currents at every scan 
step. There is a small non-reproducibility in the results of two scans, as shown in Fig. ^. The non- 
reproducibility may be attributed to a mismatch between the actual beam positions and the nominal 
ones. Therefore, it is important to check the Av and Ay values as predicted by the magnet cun^ents, 
and in particular their scales which enter linearly in the cross-section computation (Eq. |5TTD . One 
distinguishes a possible differential length scale mismatch between the two beams from a mismatch 
of their average position calibration. 

A dedicated mini-scan was performed in October where the two beams were moved in five 
equidistant steps both in x and y keeping the nominal separation between the beams constant. 
During the scan along x the beam separation was 80 |xm in x and ]xm in y. Here 80 |xm is 
approximately the width of the luminosity profile of the VDM scan (see Table This separation 
was chosen to maximize the derivative dL/dA{x), i.e. the sensitivity of the luminosity to a possible 
difference in the length scales for the two beams. If e.g. the first beam moves slightly faster than 
the second one compared to the nominal movement, the separation A(;c) gets smaller and the effect 
can be visible as an increase of the luminosity. Similarly, the beam sepai^ation used in the y scan 
was (J.m and 80 |a.m in x and y, respectively. 

The behaviour of the measured luminosity during the length-scale calibration scans is shown 
in Fig. [7|. As one can see, the points show a significant deviation from a constant. This effect may 
be attributed to different length scales of the two beams. More specifically, we assume that the real 
positions of the beams x\ 2 could be obtained from the values ;cj 2 derived from the LHC magnet 
currents by applying a correction parametrized by Ex 



and similarly for 3^1 2. The + (— ) sign in front of Ex holds for beam 1 (beam 2). Assuming a 
Gaussian shape of the luminosity dependence on A.x during the VDM scan, we get 



Here A^ = 80 \xm is the fixed nominal beam separation. A similar equation holds for the y coordi- 
nate. In the approximation of a single Gaussian shape of the beams, the width of the VDM profile, 
is defined as 



xi.2 = {y±ex/2)x\2, 



(5.3) 



1 dL 

L d{x\ +X2)/2 




(5.4) 




(5.5) 
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Figure 7. Average number of interactions (Mvelo) versus the centre of the luminous region summed over 
the twelve colliding bunches and measured during the length scale scans in x (left) and in y (right) taken 
in October. The points are indicated with small horizontal bars, the statistical errors are smaller than the 
symbol size. The straight-line fit is overlaid. 

where a^^ is width of the luminous region in the z. coordinate and Gbx is the width of beam b 
(b = 1,2). A similar equation can be written for Ly. From the slopes observed in Fig. [7| we obtain 
Ex = 2.4% and £y = —1.9%. The same behaviour is observed for all bunches separately. 

Since = (aJj — jCj) + £ (jCj +;c2)/2, the correction depends on the nominal mid-point 
between the beams {x^+X2) /2. In the first scan this nominal point was always kept at zero, there- 
fore, no correction is needed. During the second scan this point moved with nominal positions 
— )• 355.9 |J.m — )• 0. Therefore, a correction to the Ax values in Fig. is required. The central 
point should be shifted to the right (left) for the x (y) scan. The left (right) side is thus stretched 
and the opposite side is shrunk. After correction the shift between the scans is reduced in y, but 
appears now in x, so that the discrepancy cannot be fully explained by a linear coiTcction alone. 
The correction which stretches or shrinks the profiles measured in the second scan influences the 
integrals of these profiles and the resulting cross-sections very little. The latter changes on average 
by only 0.1%, which we take as an uncertainty and which we include into the systematic eiTor. In 
Table ^ the numbers are given with the correction applied. 

During a simultaneous parallel translation of both beams, the centre of the luminous region 
should follow the beam positions regardless of the bunch shapes. Since it is approximately at 
{xi +X2) /2 = [x^ +Xj)/2 and similarly for y, the corrections to the position of the centre due to 
are negligible. The luminous centre can be determined using vertices measured with the VELO. 
This provides a precise cross check of the common beam length scales (jCj +X2)/2 and {y^ +^2) /2- 
The result is shown in Fig. |S[ The LHC and VELO length scales agree within (—0.97 it 0.17)% 
and (—0.33 it 0.15)% in x and y, respectively. The scale of the transverse position measurement 
with the VELO is expected to be very precise owing to the fact that it is determined by the strip 
positions of the silicon sensors with a well-known geometry. For the cross-section determina- 
tion we took the more precise VELO length scale and multiplied the values from Table ^ by 
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Figure 8. Centre of the luminous region reconstructed with VELO tracks versus the position predicted by the 
LHC magnet currents. The points are indicated with small horizontal bars, the statistical errors are smaller 
than the symbol size. The points are fitted to a linear function. The slope calibrates the common length 
scale. 

(1 —0.0097) X (1 —0.0033) = 0.9870. In addition, we conservatively assigned a 1% systematic 
error due to the common scale uncertainty. 

In April, no dedicated length scale calibration was performed. However, a cross check is 
available from the distance between the centre of the luminous region measured with the VELO 
and the nominal centre position. The comparison of these distances between the first and second 
scan when either both beams moved symmetrically or only the first beam moved, provides a cross 
check which does not depend on the bunch shapes. From this observation the differences of the 
length scales between the nominal beam movements and the VELO reference are found to be 
(—1.3 lb 0.9)% and (1.5 it 0.9)% for and Ay, respectively. Conservatively a 2% systematic error 
is assigned to the length scale calibration for the scans taken in April. 

5.3.3 Coupling between the x and y coordinates in the LHC beams 

The LHC ring is tilted with respect to the horizontal plane, while the VELO detector is aligned 
with respect to a coordinate system where the x axis is in a horizontal plane [23]. The van der 
Meer equation (Eq. |5TTD is valid only if the particle distributions in x and in y are independent. To 
check this condition the movement of the centre of the luminous region along y is measured during 
the length scale scan in x and vice versa. This movement is compatible with the expected tilt of 
the LHC ring of 13 mrad at LHCb [23] with respect to the vertical and the horizontal axes of the 
VELO. The corresponding correction to the cross-section is negligible (< 10^^). 

To measure a possible x-y correlation in the machine the two-dimensional vertex map is studied 
by determining the centre position in one coordinate for different values of the other coordinate. 
For the analysis, data were collected with the beams colliding head-on at LHCb in Fill 1422, during 
which also the VDM scan data were taken. Figure shows the x-y profile of the luminous region. 
The centre positions of the y coordinate lie on a straight line with a slope of 79 mrad. The slopes 
found in the corresponding x-z and y-z profiles are —92 (^rad and 44 |xrad. These slopes are due 
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Figure 9. Contours of the distribution of the x-y coordinates of the luminous region. The contour Hnes show 
the values at multiples of 10% of the maximum. The points represent the >'-coordinates of the centre of the 
luminous region in different x slices. They are fitted with a linear function. 



to the known fact that the middle line between the two LHC beams is inclined with respect to the 
z axis. This is observed with beam gas events, the inclination varies slightly from fill to fill. The 
measurement of the beam directions will be described in detail in Sect. ^ Taking into account 
these known correlations of x and y with z and also the known 13 mrad tilt of the LHC ring, one 
can calculate the residual slope of the x-y correlation, which is predicted to be 77 mrad. 

If the beam profiles are two-dimensional Gaussian functions with a non-zero coiTclation be- 
tween the X and y coordinates, the cross-section relation (Eq. pTTD should be corrected. We assume 
that the x-y correlation coefficients of the two beams, ^, are similar and, therefore, close to the mea- 
sured correlation in the distribution of the vertex coordinates of the luminous region, = 0.077. 
In this case the correction to the cross-section is ^^/2 = 0.3%. We do not apply a corresponding 
correction, but instead include 0.3% uncertainty as an additional systematic error. 

5.3.4 Cross check with the z position of the luminous region 

A cross check of the width of the luminosity profile as a function of A^; is made by measuring the 
movement of the z position of the centre of the luminous region during the first VDM scan in the x 
coordinate in October (see Fig. |T^. Assuming Gaussian bunch density distributions and identical 
widths of the two colliding beams, the slope is equal to 

<iZ0 _ sin2« (^z~^x ^5 

d(Av) 4 cos^ a + a? sin^ a ' 

where Ox,; are the beam widths in the corresponding directions and dz^ is the induced shift in 
the z-coordinate of the centre of the luminous region. We approximate a, as v^ag,^. Using the 
slope observed in Fig. |T^ one gets the expected width of the luminosity profile versus A^, a/™ = 
V^Ox = 78 |a.m, in agreement with the measured value of 80 |J.m. A further cross check is described 
in Sect. \J\. 



-24- 



E 

E, 

N 

60 
40 
20 


-20 
-40 
-60 

-300 -200 -100 100 200 300 

Ax ()im) 

Figure 10. Movement of the centre of the luminous region in z during the first scan in x taken in October. 

Table 5. Cross-section of pp interactions producing at least two VELO tracks, measured in the two van der 
Meer scans in April and in October 

a^jj, (mb ) Relative uncertainty ( % ) 

April 59J 7^ 

October 58.35 3.64 

5.4 Results of the van der Meer scans 

The cross-section results obtained with the VDM scans are given in Table 0. We recall that the 
effective cross-section defined by interactions with at least two VELO tracks is used in the analysis. 
The results of the two measurement periods are consistent. During the second scan in October the 
beam movements were not continuous, so the results might suffer from hysteresis effects. In the 
second scan in April only one beam moved limiting the separation range. Therefore, only the 
first scans are shown in Table § for April and October. Both in the April and in the October 
measurements the difference observed in the results of the two scans is included as a systematic 
eiTor. The complete list of errors taken into account is given in Table ^. The uncertainties ai^e 
uncorrelated and therefore added in quadrature. As already anticipated in Sect. ^TT] the first scan 
taken in October is retained as the final result of the VDM method since it has much reduced 
systematic errors compared to the April scans. 

6. The beam-gas imaging (BGI) method 

Tracks measured in the LHCb detector allow vertices from the interactions of beam particles with 
the residual gas in the machine (beam-gas interactions) and from the beam-beam collisions to be 
reconstructed. The beam-gas imaging method is based on the measurement of the distributions 
of these vertices to obtain an image of the transverse bunch profile along the beam trajectory. 
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Table 6. Summary of relative cross-section uncertainties for the van der Meer scans in October and April. 
Due to the lower precision in the April data some systematic errors could not be evaluated and are indicated 
with 

Source relative uncertainty (%) 





October 


April 


Relative normalization stability 


0.5 


0.5 


DCCT scale 


2.7 


2.7 


DCCT baseline noise 


negligible 


3.9 


FBCT systematics 


0.2 


2.9 


Ghost charge 


0.2 


0.1 


Statistical error 


0.1 


0.9 


N\ N2 drop 


negligible 


negligible 


Emittance growth 


negligible 


negligible 


Reproducibility at nominal position 


0.4 




Difference between two scans 


2.1 


4.4 


Average beam length scale 


1.0 


2.0 


Beam length scale difference 


0.1 




Coupling of X and y coordinates 


0.3 




Total 


3.6 


7.5 



This allows the beam angles, profiles and relative positions to be determined. The residual gas in 
the beam vacuum pipe consists mainly of relatively light elements such as hydrogen, carbon and 
oxygen. 

An important prerequisite for the proper reconstruction of the bunch profiles is the transverse 
homogeneity of the visualizing gas (see Sect. |6.5.40 . A dedicated test performed in October 2010 
measured the beam-gas interaction rates as function of beam displacement in a plane peipendic- 
ular to the beam axis. The correction to the beam overlap integral due to a possible non-uniform 
transverse distribution of the residual gas is found to be smaller than 0.05% and is neglected. 

Compared to the VDM method, the disadvantage of a small rate is balanced by the advantage 
that the method is non-disruptive, without beam movements. This means that possible beam- 
beam effects are constant and effects which depend on the beam displacement, like hysteresis, are 
avoided. Furthermore, the beam-gas imaging method is applicable during physics fills. 

The half crossing-angle a is small enough to justify setting cos^ a = 1 in Eq. |0| In the ap- 
proximation of a vanishing correlation between the transverse coordinates the x and y projections 
can be factorized. At the level of precision required, the bunch shapes are well described by Gaus- 
sian distributions. Thus, their shapes are characterized in the x-y plane at the time of crossing by 
their widths a^,-, their mean positions (i = x,y), and by their bunch length Obz- The index b takes 
the values 1 and 2 according to the two beams. With these approximations, Eq. ( |1 . ID for a single 
pair of colliding bunches reduces to [24] 



NlN2f T-r 1 f 



1kJ\+ tan2 a{al + al^/{ol + a|J / af. + a|. V 2 o,^ + a, 



^2 
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where the denominator of the first factor in the product corrects for the crossing angle. The analysis 
is applied for each individual colliding bunch pair, i.e. bunch populations, event rates and beam 
profiles are considered per bunch pair. Thus, each colliding bunch pair provides an internally 
consistent measurement of the same visible cross-section. The observables a/,, and are extracted 
from the transverse distributions of the beam-gas vertices reconstructed in the bb crossings of the 
colliding bunch pairs (see Sect. 

The beam overlap-integral is then calculated from the two individual bunch profiles. The 
simultaneous imaging of the pp luminous region further constrains the beam parameters. The dis- 
tribution of /jp-coUision vertices, produced by the colliding bunch pair and identified by requiring 
— 150 < z < 150 mm, is used to measure the parameters of the luminous region. Its positions ^^-^ 
and transverse widths a^-^, 



constrain the bunch observables. Owing to the higher statistics of pp interactions compared to 
beam-gas interactions, the constraints of Eq. ( |6.2D provide the most significant input to the overlap 
integral. Equation ( |6.2D is valid only for a zero crossing angle. It will be shown in Sect. |6.5.5| that 
the approximation is justified for this analysis. 

The bunch lengths are extracted from the longitudinal distribution of the /^/^-collision ver- 
tices.^ Because the sizes Ot^ are approximately 1000 times larger than a^, the crossing angle 
reduces the luminosity by a non-negligible factor equal to the first square root factor in Eq. (|6.1D. 
The case of non-coUinear beams is described in more detail in Sect. 16.5.51 . 

The BGI method requires a vertex resolution comparable to or smaller than the transverse 
beam sizes. The knowledge of the vertex resolution is necessary to unfold the resolution from the 
measured beam profiles. The uncertainty in the resolution also plays an essential role in determin- 
ing the systematic eiTor. 

The beam-gas interaction rate determines the time needed to take a snapshot of the beam pro- 
files and the associated statistical uncertainty. When the time required to collect enough statistics 
is large compared to the time during which the beam stays stable, it becomes necessary to make 
additional corrections. This introduces systematic effects. 

6.1 Data-taking conditions 

The data used for the results described in the BGI analysis were taken in May 2010. In the data 
taken after this time the event rate was too high to select beam-gas events at the trigger level. In Oc- 
tober a more selective trigger was in place and sufficient data could be collected. However, in this 
period difficulties were observed with the DCCT data for LHC filling schemes using bunch trains. 
One should observe that the VDM data taken in October used a dedicated fill with individually 
injected bunches so that these problems were not present. 

In the selected fills, there were 2 to 13 bunches per beam in the machine. The number of 
colliding pairs at LHCb varied between 1 and 8. The trigger included a dedicated selection for 
events containing beam-gas interactions (see Sect. 0). 

^In fact, only the combination {a^, + ct^,) can be obtained. 




and a|i = 




.2 ' 



(6.2) 
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Table 7. LHC fills used in the BGI analysis. The third and fourth columns show the total number of 
(colliding) bunches ntot (wcoii)^ the fifth the typical number of particles per bunch, the sixth the period of time 
used for the analysis, and for the fills used in the BGI analysis the seventh and eighth the measured angles in 
X (in mrad) of the individual beams with respect to the LHCb reference frame (the uncertainties in the angles 
range from 1 to 5 i^xrad). The last two columns give the typical number of events per bunch used in the BGI 
vertex fits for each of the two beams. 



Fill 


part 


ntot 


Wcoll 




time (h) 


O^beam 1 




analysis 


events 1 


events 2 


1089 




2 


1 


2 X 10'" 


15 


0.209 


-0.371 


BGI 1 


1270 


720 


1090 




2 


1 


2 X 10'" 


4 


0.215 


-0.355 


BGI 2 


400 


300 


1101 




4 


2 


2 X 10'" 


6 


-0.329 


0.189 


BGI 3 


730 


400 


1104 


A 


6 


3 


2x 10'" 


5 


0.211 


-0.364 


BGI 4 


510 


350 


1104 


B 


6 


3 


2x 10'" 


5 


0.211 


-0.364 


BGI 5 


520 


350 


1117 




6 


3 


2 X 10'" 


6 


-0.327 


0.185 


BGI 6 


700 


500 


1118 




6 


3 


2x 10'" 


5 


-0.332 


0.181 


BGI 7 


500 


400 


1122 




13 


8 


2 X 10'" 


3 


-0.329 


0.182 


BGI 8 


300 


250 



The HLT runs a number of algorithms designed to select beam-gas interactions with high 
efficiency and low background. The same vertex algorithm is used for the be, eb and bb crossings, 
but different z-selection cuts are applied. For bb crossings the region —0.35 < z < 0.25 m is 
excluded to reject the overwhelming amount of interactions.^ 

6.2 Analysis and data selection procedure 

The standard vertex reconstruction algorithms in LHCb are optimized to find pp interaction vertices 
close to z = 0. This preference is removed for this particular analysis such that no explicit bias is 
present in the track and vertex selection as a function of z- The resolution of the vertex measurement 
has to be known with high precision. Details of the resolution study are given in Sect. |Ot 

The BGI method relies on the unambiguous selection of beam-gas interactions, also during bb 
crossings where an overwhelming majority of pp collisions is present. The beam-gas fraction can 
be as low as 10^^ depending on the beam conditions. The criteria to distinguish beam-gas vertices 
from pp interactions exploit the small longitudinal size of the beam spot (luminous region). As 
an additional requirement only vertices formed with exclusively forward (backward) tracks are 
accepted as beam l(2)-gas interactions and vertices are required to be made with more than ten 
tracks. A further selection on the transverse distance from the measured beam-axis is applied to 
reject spurious vertices^ (±2 mm). Due to the worsening of the resolutions for large distances from 
z = and due to the presence of pp interactions near z = 0, to determine the width of the beams 
the analysis regions are limited to —700 < z < —250 mm for beam 1 and 250 < z < 800 mm for 
beam 2.^ The selection of pp events requires —150 < z < 150 mm and only accepts vertices with 
more than 20 tracks. The background of beam-gas interactions in the pp interaction sample is 
negligible owing to the high pp event rate. 

^Due to the asymmetric VELO geometry, the background from pp vertices near the upstream end of the VELO is 
more difficult to reject, hence the asymmetric z selection. 
^Interactions in material and random associations of tracks. 

^The vertex resolution for beam 2 has a weaker z dependence, so the sensitivity is improved by enlarging the region. 
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The transverse profiles of tfie two beams are measured for each individual colliding bunch by 
projecting the vertex position on a plane perpendicular to the beam direction. The direction of the 
beam is determined on a fill-by-fill basis using the beam-gas interactions observed in be and eb 
crossings, which are free of pp interactions. The direction of the beam axis can be determined with 
1 to 5 [xrad precision depending on the fiU. 

Out of the many LHC fills only seven are selected for the BGI analysis. It is required that all 
necessary data (DCCT, FBCT, luminosity counters and vertex measurements) are present during a 
sufficiently long period and that the bunch populations and emittances are sufficiently stable during 
the selected period. The list of used fills is given in Table |7[ The table shows the total number of 
bunches and the number of bunches colliding at LHCb, the typical number of protons per bunch, 
the measured beam slopes with respect to the LHCb reference frame, and the duration of the period 
used for the analysis. 

The bunch population and size cannot be assumed to be constant during the analysis period. 
Therefore, the DCCT and FBCT data and the vertex measurements using pp interactions are binned 
in periods of 900 seconds. This binning choice is not critical. The chosen value maintains sufficient 
statistical precision while remaining sensitive to variations of the beams. The distributions of beam- 
gas interactions do not have sufficient statistics and are accumulated over the full periods shown 
in Table [7[ The analysis proceeds by determining a time-weighted average for the bunch-pair 
population product and the width and position of the pp beam spot. The weighting procedure 
solves the difficulty introduced by short periods of missing data by a logarithmic interpolation for 
the bunch populations and a linear interpolation for the bunch profiles. The averages defined by the 
latter procedure can be directly compared to the single measurement of the profiles of the single 
beams accumulated over the full period of multiple hours. A systematic eiTor is assigned to account 
for the approximations introduced by this averaging procedure. 

6.3 Vertex resolution 

The measured vertex distribution is a convolution of the true width of each beam with the resolution 
function of the detector. Since the resolution is comparable to the actual beam size (approximately 
35 (J.m in the selected fills), it is crucial to understand this effect, and to be able to unfold it from 
the reconstructed values. 

The vertex resolution is parametrized as a function of the multiplicity, or number of tracks used 
to reconstruct the vertex, and as a function of the z position of the interaction. Beam-gas vertices 
alone ai^e used to measure the positions and spatial extent of each beam; however, these events are 
rare in comparison to beam-beam vertices. To avoid binning the beam-gas vertices in both number 
of tracks and z position, beam-beam events are initially used to measure the dependence of the 
resolution on the number of tracks. Once this dependence is known, the beam-gas vertices are used 
to find the z dependence of the resolution, taking into account the contribution to the resolution 
given by the number of tracks found in each vertex. 

The resolution as a function of the number of tracks in a vertex is determined using pp in- 
teractions which occur around z = 0. The reconstructed tracks from each event are randomly split 
into two independent sets of equal size. The vertex reconstruction is applied to each set of tracks, 
and if exactly one vertex is found from each track collection it is assumed to be from the same 
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Figure 11. Primary vertex resolution ffres in the transverse directions x (full circles) and y (open circles) for 
beam-beam interactions as a function of the number of tracks in the vertex, A^rr- The curves are explained in 
the text. 



Table 8. Fit parameters for the resolution of the transverse positions x and y of reconstructed beam-beam 
interactions as a function of the number of tracks. The errors in the fit parameters are correlated. 

X y 

Factor A (mm) 0.215 ±0.020 0.202 ±0.018 

Power B 1 .023 ± 0.054 1 .008 ± 0.053 

Constant C (10^^ mm) 5.463 ± 0.675 4.875 ± 0.645 

original interaction. Then, if the number of tracks making each of these two vertices is the same, 
the residuals in x and y are taken as an estimate of the vertex resolution. 

The resolution is calculated as the width of the Gaussian function fitted to the residual dis- 
tributions divided by \/2, as there are two resolution contributions in each residual measurement. 
The resolutions are shown as a function of the number of tracks in Fig. [11]. Since the VELO is not 
fully symmetric in the x and y coordinates, the analysis is performed in both coordinates separately. 
Indeed, one observes a small, but significant difference in the resolution. The points of Fig. [TT] are 
fitted with a function which parametrizes the resolution in terms of a factor A, a power B and a 
constant C, as a function of the track multiplicity A'xr 

ares = ^±C. (6.3) 

The values of the fit parameters are given in Table |8[ 

Beam-gas vertices, selected in be and eb bunch crossings, are reconstructed in the same man- 
ner as the pp vertices. Every beam-gas event which yields two vertices after splitting the track 
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Figure 12. The z dependence of the resolution correction factor F, in x (full circles) and y (open circles) for 
beam 1-gas (left) and beam 2-gas (right) interactions. 



sample into two is used in the analysis, without requiring that the two vertices are reconstructed 
with an equal number of tracks. Accounting for the resolution coming from the track multiplicity 
means that a correction factor K. is calculated as a function of z position. This is the factor by which 
the beam-beam resolution at z = must be multiplied to find the true resolution for an event with 
a particular number of tracks, at a certain position in z- The contribution to the resolution from the 
track multiplicity is taken into account for coordinate v (v = x, j) according to 

Vl — Vl 

where the index 1,2 signifies the two vertices with measured position vi and V2, and A/^pri jA'ti^ the 
number of tracks in each. The quantities Oj^^^^ ^ are the resolutions ares expected for vertices at 
z = made of A^Tr tracks as defined in Eq. |53t 

The correction factors are plotted in Fig. |T^ In order to better understand the resolution as 
a function of z, it is instructive to consider the geometry of the VELO, shown in Fig. |T]. Figure |T^ 
shows that around the interaction point of z = the coiTcction factor is close to one, which signifies 
that the resolution is nearly independent of the type of event, whether beam-beam or beam-gas. 
This is expected, since the VELO is optimized to reconstruct vertices near z = 0. The correction 
factor increases as the vertices move away from the interaction region. 

6.4 Measurement of the beam profiles using the BGI method 

In Fig. [13| the positions of the vertices of beam-gas interactions of the single beams in be and eb 
crossings are shown in the x-z and y-z planes. The straight line fits provide the beam angles in 
the corresponding planes. Whereas we can use the non-colliding bunches to determine the beam 
directions, the colliding bunches are the only relevant ones for luminosity measurements. 

As an example the x and y profiles of one colliding bunch pair are shown in Fig. [n[ The 
physical bunch size is obtained after deconvolving the vertex resolution. The resolution function 
and physical beam profile are drawn separately to show the importance of the knowledge of the 
resolution. In Fig. [T5| the corresponding fits to the luminous region of the same bunch pair are 
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Figure 13. Positions of reconstructed beam-gas interaction vertices for be (black points) and eb (grey points) 
crossings during Fill 1104. The measured beam angles ai,2 and offsets 5i.2 at z = in the horizontal (top) 
and vertical (bottom) planes are shown in the figure. 



shown, both for the full fill duration and for a short period of 900 s. The fits to the distributions 
in X and y of the full fill have a ;^^/ndf 10, probably due to the emittance growth of the beam. 
The corresponding fits to data taken during the shorter period of 900 s give satisfactory values, 
;^^/ndf 1. The fact that the resolution at z = is small compared to the size of the luminous 
region makes it possible to reach small systematic uncertainties in the luminosity determination as 
shown in Sect. 1531 

For non-colliding bunches it is possible to measure the width of the beam in the region of the 
interaction point at z = since there is no background from pp collisions. One can compare the 
measurement at the IP with the measurement in the region outside the IP which needs to be used 
for the colliding bunches. After correcting for the resolution no difference is observed, as expected 
from the values of j8* of the beam optics used during the data taking. One can also compare 
the width measurements of the colliding bunches far from the IP using the beam-gas events from 
beam 1 and beam 2 to predict the width of the luminous region using Eq. |5^ Figure |T^ shows 
that there is overall consistency. In addition to the data used in the BGI analysis described here, 
also higher statistics data from later fills are used for this comparison. The cross check reaches 
a precision of 1-1.6% for the consistency of the width measurements at large z compared to the 
measurement at z = 0, providing good evidence for the correctness of the parametrization of the z 
dependence of the vertex resolution. 

The relations of Eq. are used to constrain the width and position measurements of the 
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Figure 14. Distributions of the vertex positions of beam-gas events for beam 1 (top) and beam 2 (bottom) 
for one single bunch pair (ID 2186) in Fill 1 104. The left (right) panel shows the distribution in x (y). The 
Gaussian fit to the measured vertex positions is shown as a solid black curve together with the resolution 
function (dashed) and the unfolded beam profile (shaded). Note the variable scale of the horizontal axis. 



single beams and of the luminous region in botli coordinates separately. Given the high statistics of 
vertices of the luminous region the pp events have the largest weight in the luminosity calculation. 
Effectively, the beam-gas measurements determine the relative offsets of the two beams and their 
width ratio p 

pi = 02il 0\i (/ = x,y). (6.5) 

According to Eq. |57T1 , neglecting crossing angle effects and beam offsets, the luminosity is 
proportional to A~j^, 

= I \ ■ (6-6) 

This quantity can be rewritten using the definition of p, and Eq. 
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Figure 15. Distributions of vert;ex positions of pp interactions for the full fill duration (left) and for a 900 s 
period in the middle of the fill (right) for one colUding bunch pair (ID 2186) in Fill 1 104. The top, middle 
and bottom panels show the distributions in x, >' and z, respectively. The Gaussian fit to the measured vertex 
positions is shown as a solid black curve together with the resolution function (dashed) and the unfolded 
luminous region (shaded). Owing to the good resolution the shaded curves are close to the soUd curves and 
are therefore not clearly visible in the figures. The fit to the z coordinate neglects the vertex resolution. Note 
the variable scale of the horizontal axis. 
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Figure 16. Comparison of the prediction for the luminous region width from measurements based on beam- 
gas events of individual bunches which are part of a colliding bunch pair with the direct measurement of the 
luminous region width for these colliding bunches. The panels on the left show the results for bunches in the 
fills with j3* = 2 m optics used in this analysis, the right panels show four colliding bunches in a fill taken 
with j3* = 3.5 m optics. The fill and bunch numbers are shown on the vertical axis. The vertical dotted line 
indicates the average and the solid lines the standard deviation of the data points. The lowest point indicates 
the weighted average of the individual measurements; its error bar represents the corresponding uncertainty 
in the average. The same information is given above the data points. The fills with the )3* = 3.5 m optics are 
not used for the analysis due to the fact that larger uncertainties in the DCCT calibration were observed. 

which shows, especially for nearly equal beam sizes, the weight of the measurement of the width 
of the luminous region cr^j in the luminosity determination. The expression has its maximum for 
Pi = 1 which minimizes the importance of the measurement errors of p,. 

The luminosity changes if there is an offset between the two colliding bunches in one of the 
transverse coordinates. An extra factor appeared already in Eq. ^7T] to take this into account 

Corset =n-P (4^1^), (6.8) 
i=x,y V ^ "iM "2/ / 

which is unity for head-on beams. 
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By examining both relations of Eq. ^3 ^ system of two constraint equations and six measur- 
able quantities emerges which can be used to improve the precision for each transverse coordinate 
separately. This fact is exploited in a combined fit where the individual beam widths On, G2i, and 
the luminous region width cr^j together with the corresponding position values and are used 
as input measurements. Several choices are possible for the set of four fit-parameters, trivially the 
set On, 02i, ^li, can be used. The set £, (Eq. |53[ ), p, (Eq. |63| ), A^,- = — and is used 
which makes it easier to evaluate the corresponding luminosity enor propagation. The results for 
the central values are identical, independently of the set used. 



6.5 Corrections and systematic errors 

In the following, corrections and systematic error sources affecting the BGI analysis will be de- 
scribed. The uncertainties related to the bunch population normalization have already been dis- 
cussed in Sect. 3. 



6.5.1 Vertex resolution 

As mentioned above, the uncertainty of the resolution potentially induces a significant systematic 
error in the luminosity measurement. To quantify its effect the fits to the beam profiles have been 
made with different choices for the parameters of the resolution functions within the limits of their 
uncertainty. One way to estimate the uncertainty is by comparing the resolution of simulated pp 
collision events determined using the method described in Sect. i.e. by dividing the tracks into 
two groups, with the true resolution which is only accessible in the simulation owing to the prior 
knowledge of the vertex position. The uncertainty in the number of tracks (A'n) dependence at 
z = is estimated in this way to be no larger than 5%. 

The uncertainty in the z dependence is estimated by analysing events in eb and be crossings. 
For these crossings all events, including the ones near z = 0, can be assumed to originate from 
beam-gas interactions. A cross check of the resolution is obtained by comparing the measurement 
of the transverse beam profiles at the z range used in bb events with the ones obtained near z = 0. 
Together with the comparison of the width of the luminous region and its prediction from the beam- 
gas events we estimate a 10% uncertainty in the z dependence of the resolution. It should be noted 
that the focusing effect of the beams towards the interaction point is negligible compai^ed to the 
precision needed to determine the resolution. 

The effect of the uncertainties in the A^Tr dependence and z dependence on the final results are 
estimated by repeating the analysis varying the resolution within its uncertainty. The conservative 
approach is taken to vary the different dependencies coherently in both x and y and for the depen- 
dence on Njr and z simultaneously. The resulting uncertainties in the cross-section depend on the 
widths of the distributions and are therefore different for each analysed bunch pair. 

6.5.2 Time dependence and stability 

The beam and data-taking stability are taken into account when selecting suitable fills to perform 
the beam-gas imaging analysis. This is an essential requirement given the long integration times 
needed to collect sufficient statistics for the beam-gas interactions. A clear decay of the bunch- 
populations and emittance growth is observed over these long time periods. It is checked that these 
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variations are smooth and that the time-average is a good approximation to be used to compare 
with the average beam profiles measured with vertices of beam-gas interactions. No significant 
movement of the luminous region is observed during the fills selected for the BGI analysis. The 
systematics introduced by these variations are minimized by the interpolation procedure described 
in Sect. and are estimated to amount to less than 1%. 

6.5.3 Bias due to unequal beam sizes and beam offsets 

When the colliding bunches in a pair have similar widths (p,- = 1), the p -dependence in Eq. |577| , 
p,7(l + pf ), is close to its maximum. Thus, when the precision of measuring p is similar to its 
difference from unity, the experimental estimate of the p-factor is biased towards smaller values. In 
the present case the deviation from unity is compatible with the statistical error of the measurement 
for each colliding bunch pair. These values are typically 15% in the x coordinate and 10% in the y 
coordinate. The size of the "p bias" effect is of the order of 1% in ;c and 0.5% in y. 

A similar situation occurs for the offset factor Coffset for bunches colliding with non-zero rela- 
tive transverse offset. The offsets are also in this case compatible with zero within their statistical 
errors and the correction can only take values smaller than one. The average expected "offset bias" 
is typically 0.5% per transverse coordinate. 

Since these four sources of bias (unequal beam sizes and offsets in both transverse coordinates) 
act in the same direction, their overall effect is no longer negligible and is corrected for on a bunch- 
by-bunch basis. We assume a systematic error equal to half of the connection, i.e. typically 1.5%. 
The correction and associated uncertainty depends on the measured central value and its statistical 
precision and therefore varies per fill. 

6.5.4 Gas pressure gradient 

The basic assumption of the BGI method is the fact that the residual gas pressure is uniform in 
the plane transverse to the beam direction and hence the interactions of the beams with the gas 
produce an image of the beam profile. An experimental verification of this assumption is performed 
by displacing the beams and recording the rate of beam-gas interactions at these different beam 
positions. In Fill 1422 the beam was displaced in the x coordinate by a maximum of 0.3 mm. 
Assuming a linear behaviour, the upper limit on the gradient of the interaction rate is 0.62 Hz/mm 
at 95% CL compared to a rate of 2.14 it 0.05 Hz observed with the beam at its nominal position. 
When the profiles of beam 1-gas and beam 2-gas interactions are used directly to determine the 
overlap integral Agff, the relative error 5a on the overlap integral is given by 

Aeff(a = 0) a^a^ 

where a is the gradient, b the rate when the beam is at its nominal position, and cJt is the true beam 
width, using the approximation of equal beam sizes. This result has been derived by comparing 
the overlap integral for beam images distorted by a linear pressure gradient with the one obtained 
with ideal beam images. With the measured limit on the gradient, the maximum relative effect on 
the overlap is then estimated to be less than 4.2 x 10^'*. However, the BGI method uses the width 
of the luminous region measured using pp interactions as a constraint. This measurement does not 
depend on the gas pressure gradient. The gas pressure gradient enters through the measurements of 
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the individual widths which are mainly used to determine the ratio between the two beam widths. 
These are equally affected, thus, the overall effect of an eventual gas pressure gradient is much 
smaller that the estimate from Eq. and can safely be neglected in the analysis. 



6.5.5 Crossing angle effects 

The expression for the luminosity (Eq. |OD contains a correction factor for the crossing angle Ca 
of the form [24] 

Ca = [1 +tan2 a{ol + ol)/{ol + a|,)] . (6.10) 

For a vanishing crossing angle and equal bunch lengths, the bunch length is obtained from the 
beam spot measurement assuming that the two beams have equal size, by a, = \fl.o^^. In the 
presence of a crossing angle the measured length of the luminous region depends on the lengths of 
the bunches, on the crossing angle and on the transverse widths of the two beams in the plane of the 
crossing angle. The bunch lengths need not necessarily be equal. Evaluating the overlap integral 
of the two colliding bunches over the duration of the bunch crossing, one finds for the width of the 
luminous region in the z coordinate 



tan^ a 4 cos^ a 



(6.11) 



Solving for af, + a2,, 
quantities a, a^^, a^^, Ou, and a2x 



the right-hand side of Eq. 6.10 can be written in terms of the measured 



Ca 



1 + - 



4sin^a(7L 



(6.12) 



(1 - (tanac7^,/cy^J2) (^^2^ + ^^2^ 

The dependence of the estimate of a^^ on a, and of the overall correction on Ox is shown in Fig. |T7| 
for typical values of the parameters. The difference with respect to a naive calculation assuming 
equal beam sizes and using the simple \/2 factor to obtain the bunch lengths from the luminous 
region length is in all relevant cases smaller than 1%. For the beam conditions in May 2010 the 
value of the crossing angle correction factor Ca is about 0.95. To take into account the accuracy of 
the calculation a 1 % systematic error is conservatively assigned to this factor. 

There are other small effects introduced by the beam angles. The average angle of the beams 
is different from in the LHCb coordinate system. This small difference introduces a broadening 
of the measured transverse widths of the luminous region, since in this case the projection is taken 
along the nominal LHCb axis. Another effect is more subtle. The expression (Eq. |0| ) for the 
width of the luminous region assumes a vanishing crossing angle. It is still valid for any crossing 
angle if one considers the width for a fixed value of z- When applying Eq. ^3 ^ function of z 
one can show that the centre of the luminous region is offset if in the presence of a non-vanishing 
crossing angle the widths of the two beams are not equal. Thus, when these two conditions are met 
the luminous region is rotated. The rotation angle (^,- (/ = x,y) is given by 



tan = tan a. 



(6.13) 



where p is defined in Eq. 
is smaller than 10^^. 



1+pr 

With the parameters observed in this analysis the effect of the rotation 
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Figure 17. Left: the dependence of the length of the luminous region a^^ on the single bunch length d, 
under the assumption that both beams have equal length bunches. The dotted line shows the \/2 behaviour 
expected in the absence of a crossing angle. The solid black line shows the dependence for equal transverse 
beam sizes (7^ = 0.045 mm, the shaded region shows the change for p = 1.2 keeping the average size 
constant. Right: the dependence of the luminosity reduction factor Ca on the transverse width of the beam 
Ox for a value of ffgj^ = 35 mm. The solid line shows the full calculation for p — I (equal beam widths) with 
the shaded area the change of the value up to p = 1 .2, keeping the transverse luminous region size constant. 
The dotted line shows the result of the naive calculation assuming a simple y/l relation for the length of the 
individual beams. All graphs are calculated for a half crossing-angle a = 0.2515 mrad. 

6.6 Results of the beam-gas imaging method 

With the use of the beam-gas imaging method seven independent measurements of an effective 
reference cross-section were performed. The main uncertainties contributing to the overall preci- 
sion of the cross-section measurement come from the overlap integral and from the measurement 
of the product of the bunch populations. The systematic error in the overlap integral is composed 
of the effect of the resolution uncertainties, the treatment of the time dependence, the treatment 
of the bias due to the non-linear dependencies in p and A((^) and the crossing angle corrections. 
It also takes into account small deviations of the beam shape from a single Gaussian observed in 
the VDM scans. The normalization error has components from the DCCT scale, and its baseline 
fluctuations, FBCT systematics, and systematics in the relative normalization procedure (Sect. 
For multi-bunch fills the results obtained for each colliding pair are first averaged taking the coitc- 
lations into account. The results of the averaging procedure applied on a per-fill basis are shown 
in Table ^ For fills with multiple bunches the numbers are a result of an average over individual 
bunches. Errors are divided into two types: correlated and uncorrelated errors. On a fiU-by-fiU 
basis the statistical errors, ghost charge and DCCT baseUne corrections are treated as uncorrelated 
errors. The latter two sources are of course correlated when bunches within one fill are combined. 
The FBCT systematic uncertainty, which is dominated by the uncertainty in its offset is treated 
taking into account the fact that the sum is constrained. Owing to the constraint on the total beam 
current provided by the DCCT, averaging results for different colliding bunch pairs within one fill 
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Table 9. Measurements of the cross-section ff^j^ with the BGI method per fill and overall average (third 
column). All errors are quoted as percent of the cross-section values. DCCT scale, DCCT baseline noise, 
FBCT systematics and Ghost charge are combined into the overall Beam normalization error. The Width syst 
row is the combination of Resolution syst (the systematic error in the vertex resolution for pp and beam-gas 
events), Time dep. syst (treatment of time-dependence) and Bias syst (unequal beam sizes and beam offset 
biases), and is combined with Crossing angle (uncertainties in the crossing angle correction) into Overlap 
syst. The Total error is the combination of Relative normalization stability. Beam normalization. Statistical 
error, and Overlap syst. Total systematics is the combination of the latter three only and can be broken down 
into Uncorrelated syst and Correlated syst, where "uncorrelated" applies to the avergaging of different fills. 
Finally, Excluding norm is the uncertainty excluding the overall DCCT scale uncertainty. The grouping 
of the systematic errors into (partial) sums is expressed as an indentation in the first column of the table. 
The error components are labelled in the second column by u, c or f dependending on whether they are 
uncorrelated, fully correlated or correlated within one fill, respectively. 
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Width syst 




3.20 


3.18 


3.43 


3.05 


3.56 


2.83 


2.99 


3.34 


Crossing angle 


c 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


Overlap syst 




3.35 


3.33 


3.58 


3.21 


3.70 


3.00 


3.15 


3.49 


Uncorrelated syst 
Correlated syst 
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f 
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4.35 
4.45 


3.08 
4.43 
5.39 


3.07 
4.62 
5.55 


2.80 
4.25 
5.08 


2.18 
4.62 
5.11 


2.44 
4.08 
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2.47 
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2.01 
4.44 
4.88 
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6.75 
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5.01 
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4.31 
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4.82 
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Figure 18. Results of the beam-gas imaging method for the visible cross-section of the pp interactions 
producing at least two VELO tracks, CJ^j^. The results for each fill (indicated on the vertical axis) are obtained 
by averaging over all colliding bunch pairs. The small vertical lines on the error bars indicate the size of the 
uncorrelated errors, while the length of the error bars show the total error The dashed vertical line indicates 
the average of the data points and the dotted vertical lines show the one standard-deviation interval. The 
weighted average is represented by the lowest data point, where the error bar corresponds to its total error. 



reduces the error introduced by the FBCT offset uncertainty. A usual error propagation is applied 
taking the inverse square of the uncorrelated errors as the weights. The difference with respect to 
the procedure applied for the VDM method is due to the fact that a fit using the FBCT offsets as 
free parameters cannot be applied here. Not all fills have the required number of crossing bunch 
pairs, and the uncorrelated bunch-to-bunch errors are too large to obtain a meaningful result for the 
FBCT offset. 

Each colliding bunch pair provides a self-consistent effective cross-section determination. The 
analysis proceeds by first averaging over all individual colliding bunch pairs within a fill and then 
by averaging over fills taking all correlations into account. Thus, an effective cross-section result 
can also be quoted per fill. These are shown in Fig. |TE|. The spread in the results is in good 
agreement with the expectations from the uncorrelated errors. 

The final beam-gas result for the effective cross-section is: a^j^ = 59.9 ± 2.7 mb. The un- 
certainties in the DCCT scale error and the systematics of the relative normalization procedure are 
in common with the VDM method. The uncertainty in a^jj, from the BGI method without these 
common errors is 2.2 mb. 
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Figure 19. Visible cross-section measurement using the beam-beam imaging method for twelve different 
bunch pairs (filled circles) compared to the cross-section measurements using the VDM method (open cir- 
cles). The horizontal line represents the average of the twelve beam-beam imaging points. The error bars 
are statistical only and neglect the correlations between the measurements of the profiles of two beams. The 
band corresponds to a one sigma variation of the vertex resolution parameters. 



7. Cross checks with the beam-beam imaging method 

During the VDM scan the transverse beam images can be reconstructed with a method described 
in Ref. [15]. When one beam (e.g. beam 1) scans across the other (beam 2), the differences 
of the measured coordinates of pp vertices with respect to the nominal position of beam 2 are 
accumulated. These differences are projected onto a plane transverse to beam 2 and summed over 
all scan points. The resulting distribution represents the density profile of beam 2 when the number 
of scan steps is large enough and the step size is small enough. By inverting the procedure, beam 1 
can be imaged using the relative movement of beam 2. Since the distributions are obtained using 
measured vertex positions, they are convolved with the corresponding vertex resolution. After 
deconvolving the vertex distributions with the transverse vertex resolution a measurement of the 
transverse beam image is obtained. This approach is complementary to the BGl and VDM methods. 

The beam-beam imaging method is applied to the first October VDM scan since the number of 
scanning steps of that scan is twice as large as that of the second scan. The events are selected by the 
minimum bias trigger. Contrary to the 22.5 kHz random trigger events with only luminosity data, 
they contain complete detector information needed for the tracking. The minimum-bias trigger 
rate was limited to about 130 Hz on average. The bias due to the rate limitation is corrected 
by normalizing the vertex distributions at every scan point to the measured average number of 
interactions per crossing. 

The transverse planes with respect to the beams are defined using the known half crossing- 
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angle of 170 |xrad and the measured inclination of the luminous ellipsoid with respect to the z axis 
as discussed in Sect. |5.3.3| . The measured common length scale correction and the difference in 
the length scales of the two beams described by the asymmetry parameters Ex^y in Eq. p3| , are also 
taken into account. 

The luminosity overlap integral is calculated numerically from the reconstructed normalized 
beam profiles. The effect of the VELO smearing is measured and subtracted by comparing with 
the case when the smearing is doubled. The extra smearing is performed on an event-by-event 
basis using the description of the resolution given in Eq. with the parameters A, B and C taken 
from Table ^ To improve the vertex resolution, only vertices made with more than 25 tracks are 
considered. This reduces the average cross-section correction due to the VELO resolution to 3.7%. 

Similar to the BGI method, the beam-beam imaging method measures the beam profiles per- 
pendicular to the beam directions. For the luminosity determination in the presence of a crossing 
angle their overlap should be corrected by the factor Ca (see Eq. |6.1(J| ) due to a contribution from 
the length of the bunches. The average correction for the conditions during the VDM fill in Octo- 
ber is 2.6%. A bunch-by-bunch comparison of the cross-section measurement with the beam-beam 
imaging method and the VDM method is shown in Fig. |T^ The cross-section is measured at the 
nominal beam positions. The FBCT bunch populations with zero offsets normalized to the DCCT 
values and corrected for the ghost charge are used for the cross-section determination. The band in- 
dicates the variation obtained by changing the vertex resolution by one standard deviation in either 
direction. The obtained cross-section of 59.1 mb is in good agreement with the value of 58.4 mb 
reported in Table ^ The comparison is very sensitive since the overall bunch population normal- 
ization and the length scale uncertainty are in common. Uncorrelated errors amount to about 1%. 
The main uncorrelated errors in the beam-beam imaging method are the VELO systematics and 
the statistical error which are each at the level of 0.4%. The main uncorrelated errors in the VDM 
method are the stability of the working point (0.4%) and the statistics (0.1%). The difference be- 
tween the two methods is smaller than but similar to the difference between the two October scan 
results observed with the VDM method. 

The width of the VDM rate profile and the widths of the individual beams are related following 
Eq. |53[ Thus the widths Lx.y are directly measured with the VDM rate profile and predicted 
using the measured widths of the beams. The widths can be compared directly using the RMS 
of the distributions, the widths (RMS) of single Gaussian fits or the RMS of double Gaussian 
fits. The variation among these different values are of the order of 1% and limit the sensitivity. 
However, it should be noted that these are just numerical differences; Eq. holds for arbitrary 
beam shapes. The ratio of the measured and predicted width is 0.994 and 0.996 in the x and y 
coordinate, respectively. The statistical uncertainties are 0.3% and the uncertainties due to the 
knowledge of the vertex resolution are 0.2%. Considering the sensitivity of the comparison, we 
note good agreement. 

8. Results and conclusions 

The beam-gas imaging method is applied to data collected by LHCb in May 2010 using the residual 
gas pressure and provides an absolute luminosity normalization with a relative uncertainty of 4.6%, 
dominated by the knowledge of the bunch populations. The measured effective cross-section is in 
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Table 10. Averaging of the VDM and BGI results and additonal uncertainties when applied to data-sets used 
in physics analyses. 





Average 


VDM 


BGI 


Cross-section (mb) 


58.8 


58.4 


59.9 


DCCT scale uncertainty (%) 


2.7 


2.7 


2.7 


Uncorrelated uncertainty (%) 


2.0 


2.4 


3.7 


Cross-section uncertainty (%) 


3.4 


3.6 


4.6 


Relative normalization stability (%) 


0.5 






Use of average value of jj.^.^^ (%) 


0.5 






Additional uncertainty for other data-sets (%) 


0.7 






Total uncertainty for large data sets (%) 


3.5 







agreement with the measurement performed with the van der Meer scan method using dedicated 
fills in April 2010 and October 2010. The VDM method has an overall relative uncertainty of 3.6%. 
The final VDM result is based on the October data alone which give significantly lower systematic 
uncertainties. The common DCCT scale error represents a large part of the overall uncertainty 
for the results of both methods and is equal to 2.7%. To determine the average of the two results 
the common scale should be removed before calculating the relative weights. Table |T^ shows the 
ingredients and results of the averaging procedure. The combined result has a 3.4% relative error. 

Since the data-sets used for physics analysis contain only a subset of all available information 
(see Sect. a small additional error is introduced e.g. by using jj.^^^ information averaged over 
bunch crossings. Together with the uncertainty introduced by the long term stability of the rela- 
tive normalization this results in a final uncertainty in the integrated luminosity determination of 
3.5%. We have taken the conservative approach to assign a 0.5% uncertainty representing the rel- 
ative normalization variation to all data-sets and not to single out one specific period as reference. 
The results of the absolute luminosity measurements are expressed as a calibration of the visible 
cross-section a^jj,. This calibration has been used to determine the inclusive cross-section in pp 
collisions at ^/s = 7 TeV [25].^ 

The relative normalization and its stability have been studied for the data taken with LHCb in 
2010 (see Sect. 0). Before the normalization can be used for other data-sets an appropriate study 
of the relative normalization stability needs to be performed. 

While the VDM data have been taken during dedicated fills, no dedicated data taking periods 
have yet been set aside for the BGI method. It is, therefore, remarkable that this method can reach 
a comparable precision. A significantly improved precision in the DCCT scale can be expected 
in the near future. In addition, a controlled pressure bump in the LHCb interaction region would 
allow us to apply the beam-gas imaging method in a shorter period, at the same time decreasing the 
effects from non-reproducibility of beam conditions and increasing the statistical precision. The 
main uncertainty in the VDM result, apart from the scale error, is due to the lack of reproducibility 
found between different scanning strategies. Dedicated tests will have to be designed to understand 

''in fact, for the early data-taking period on which this measurement is based, the hit count in the SPD is used to 
define the visible cross-section. This cross-section differs from CT^j^ defined in this paper by 0.5%. 
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these differences better. Finally, it is also very advantageous to perform beam-gas measurements in 
the same fill as the van der Meer scans. This would allow cross checks to be made with a precision 
which does not suffer from scale uncertainties in the bunch population measurement. Furthermore, 
a number of parameters which Umit the precision of the BGI method can be constrained indepen- 
dently using the VDM scan data, such as the relative beam positions. To improve the result of the 
BGI method a relatively large jS* value should be chosen, such as 10 m. The precision of the VDM 
method does not, in principle, depend on j3*. 
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